Decoherence of a quantum memory coupled to a collective spin bath 
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We study the quantum dynamics of a single qubit coupled to a bath of interacting spins as a 
model for decoherence in solid state quantum memories. The spin bath is described by the Lipkin- 
Meshkov-Glick model and the bath spins are subjected to a transverse magnetic field. We investigate 
the qubit interacting via either an Ising- or an XY-type coupling term to subsets of bath spins of 
differing size. The large degree of symmetry of the bath allows us to find parameter regimes where 
the initial qubit state is revived at well defined times after the qubit preparation. These times may 
become independent of the bath size for large baths and thus enable faithful qubit storage even 
in the presence of strong coupling to a bath. We analyze a large range of parameters and identify 
those which are best suited for quantum memories. In general we find that a small number of 
links between qubit and bath spins leads to less decoherence and that systems with Ising coupling 
between qubit and bath spins are preferable. 
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I. INTRODUCTION 

The transfer and storage of quantum information be- 
tween different physical systems is of crucial importance 
to quantum communication and distributed quantum 
computing I]. Experimentally quantum communication 
is one of the most advanced areas of quantum technology 
with promising commercial applications being a realistic 
possibility in the near future [2J. However, for its full po- 
tential to be realized over relevant distances, limitations 
caused by noisy transmission lines, e.g. the exponential 
scaling of photon losses in an optical fiber with its length, 
need to be overcome. 

Several promising proposals to solve this problem ex- 
ist. The most well studied is the use of quantum re- 
peaters which segment the transmission line into shorter 
pieces 0,13]. By applying sophisticated entanglement 
purification schemes and entanglement swapping at these 
repeater units, high-fidelity entangled pairs can be estab- 
lished over much larger distances than direct transmis- 
sion could feasibly permit. A recent alternative method 
has been proposed that is based on the idea of entan- 
glement percolation in a quantum network. Given nodes 
that are initially connected by partially entangled pure 
states @ or mixed states 0, schemes based on classical 
bond percolation have been shown to enable the creation 
of maximally entangled singlet states between arbitrary 
points in the network with aprobability independent of 
the distance between them [7]. Both quantum repeater 
and percolation schemes require that quantum informa- 
tion is stored locally at the repeaters or nodes of the net- 
work and crucially their operation degrades in the pres- 
ence of decoherence in these memories. 
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In this work we investigate the quantum dynamics of 
a single qubit quantum memory coupled to an environ- 
ment of interacting spin-1/2 particles |8j ■ Such a model of 
decoherence has relevance for solid-state quantum mem- 
ories involving the nuclear spin [9|, llCf , electron spin in a 
semiconducting quantum dot and nitrogen vacancy 
centers in diamond [l2j . Previous studies examining the 
decoherence induced by spin baths have mostly consid- 
ered the so-called 'central spin model' in which the qubit 
is coupled isotropically to all spins in the bath. Early 
work has analyzed the decoherence due to independent 
bath spins whilst most studies since have inves- 

tigated ID models with nearest-neighbor couplings [l7r - 
|22| ]. The effect of an infinite-range interaction amongst 
bath spins has also been studied within the central spin 
model I23W26II . However, as pointed out by Rossini et al. 
in Ref. [271 ]. the assumption of a central spin may not 
be valid for many physical systems. Thus, in their work 
j27| the authors depart from the central spin model and 
consider just a few links between the qubit and a ID spin 
chain. Here we consider the spin bath as being described 
by the Lipkin-Meshkov-Glick (LMG) model (see Eq. |3J}, 
which possesses infinite-range interactions, and similarly 
go beyond the central spin model by examining the deco- 
herence of a single qubit as its exposure to the spin bath 
is varied from coupling to just one bath spin through to 
interacting with all bath spins. 

The LMG spin bath represents an ideal benchmark for 
investigating the effects of the varying qubit interactions 
with a non-trivial spin bath. Firstly, it possesses per- 
mutational symmetry, which can be exploited to exactly 
solve its dynamics numerically for very large number of 
spins (~ 100-1000). This in itself is a rare property for 
many-body quantum systems. The large degree of sym- 
metry in the bath also allows us to identify parameter 
regimes where highly non-Markovian [15} coherence prop- 
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erties are displayed within the quantum memory. We 
quantify this by its rephasing time and periodic entan- 
glement breaking, which are found to be independent of 
the bath size. Secondly, the LMG model has also at- 
tracted much attention [28-32] for exhibiting either first- 
or second-order equilibrium quantum phase transitions 
dependent on the type of intra-bath coupling. As is now 
well known the entanglement properties of the ground 
state are strongly affected at criticality [33Tj36j and this 
in turn can have a profound effect on the induced deco- 
herence of a coupled qubit [TtJ • We explore the effects of 
criticality and in a similar fashion to H. T. Quan et al. 
in Ref. [17| our findings point to a possible use of a cou- 
pled quantum memory as an apparatus for the dynamical 
measurement of a QPT through its coherence properties. 

The structure of this paper is as follows. In Sec. |TT] 
we introduce the model and the bath Hamiltonian, and 
describe the ground state properties of the bath under 
various parameter regimes. We also introduce the de- 
phasing Ising interaction and dissipative XY interaction 
terms and discuss the measures used to quantify deco- 
herence in our system. The results are presented in Sec. 
IIIII and IIVI for dephasing and dissipative interactions re- 
spectively, with conclusions drawn in Sec. [Vj 



II. THE MODEL 

The model under consideration is that of a single spin- 
1/2 (qubit) Q interacting with a spin bath R. The pur- 
pose of the investigation is two-fold: the primary focus 
is to study how the interaction with the bath affects the 
evolution of Q, however in doing so we also use Q as a 
probe to study the bath dynamics across criticality. The 
Hamiltonian of the global system, depicted schematically 
in Fig. [T] is of the standard form 



H, 



H r + Ht 



(1) 



where Hq and Hr are the free Hamiltonians of the qubit 
and the bath respectively, and Hi is the interaction term. 
For the qubit we assume 



H, 



Q 



(2) 



where q z is the Pauli z-operator and uj is the energy dif- 
ference between the ground and excited states; we work 
in units with h = 1. In the following subsections, we 
will introduce the bath Hamiltonian Hr — Hlmg an d 
the two different interaction terms under consideration. 
Throughout this work, the subsystems Q and R will be 
initiated at time t = in a product state, with the bath 
in its ground state. 



A. The Lipkin-Meshkov-Glick bath 

The bath is described by the Lipkin-Meshkov-Glick 
model, which was originally proposed in nuclear physics 




FIG. 1. (Color on line) A schematic of the model under inves- 
tigation. A single spin-1/2 is coupled via k links (red lines) 
to a completely connected spin bath R. Note that the lengths 
of the lines do not correspond to the strength of the interac- 
tions. The green dash-dotted line indicates a particular bath 
partition. 



to describe shape phase transitions in nuclei |37|. We 
note that it has since found relevance in describing the 
Josephson effect in two-mode Bose-Einstein condensates 
38] and a dissipative implementation of it has also re- 
cently been proposed in optical cavity QED [39]. The 
model involves a system of N spins- 1/2 mutually inter- 
acting in the x-y plane and embedded in a transverse 
magnetic field. The Hamiltonian is given by 



A 



N 



H LMG = £ (ofo? + T^jO - h^af 



N 

t<j 

-^(S 2 + 7 S 2 )-2hS z + 



-(1 + 7), (3) 



where S a — J2i=i fJ ? is the collective spin operator 
for the direction a = x,y,z; and the a a are the Pauli 
operators for spins in the bath. The quantities A, h, and 7 
characterize the interaction strength between bath spins, 
the magnetic field strength, and the anisotropy in the x- 
y plane, respectively. The above Hamiltonian describes 
an infinite-ranged XY model, with the prefactor 1 /N on 
the first term necessary to ensure a finite free energy per 
spin in the thermodynamical limit. The cases 7=1 and 
7 = correspond to infinite-ranged XX and transverse 
Ising models respectively. In the following, we consider 
only ferromagnetic interactions A > within the bath for 
anisotropy values < 7 < 1. 

The mutual coupling of all spins described by Hlmg 
results in the energy eigenstates of the free bath being 
invariant under particle exchange. Furthermore, Hlmg 
commutes with the total spin operator S 2 = 5 2 + 5 2 + 5 2 , 
and so the energy eigenstates are divided into symmet- 
ric subspaces Ss defined by the spin quantum number 5 
of S 2 and each of dimensionality 25+1. This symme- 
try vastly reduces the complexity of studies of the ground 
state (GS) of the bath. For ferromagnetic intra-bath cou- 
plings the GS lies in the maximal spin subspace §jv/2i and 
so its calculation is an eigenvalue problem scaling as just 
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FIG. 2. (Color on line) Upper: the energy spectrum for 
isotropic (a) and anisotropic (b) ferromagnetically coupled 
baths of N — 30 spins. Blue lines indicate levels in the sub- 
space §jv/2, green lines the subspace §nv/2)-li an d red lines 
the subspace §(jv72)-2- Lower: GS magnetization (per spin) 
in the z-direction for isotropic (c) and anisotropic (d) ferro- 
magnetically coupled baths of N = 500 spins. 



2V + 1. 

For isotropic intra-bath coupling (7=1) Hlmg also 
commutes with S z and is thus diagonal in the symmetric 
eigenbasis spanned by Dicke states \S,M) (M is the 
quantum number for S z ). Dicke states are degenerate 
except for those belonging to the maximal spin subspace 
§/v/2- For an ensemble of N spins, Dicke states in the 
subspace Sjv/2 are entangled states except for the two 
spin-polarized states \N/2,N/2) and \N/2, -N/2). 



1. Ground state properties 

The ground state of the ferromagnetic bath lies in the 
maximal spin subspace §jv/2 for all 7 and A > 0. A mean- 
field approach [4ll |42| predicts a second-order quantum 
phase transition in the thermodynamical limit at the crit- 
ical point A = h for all 7 > 0. 

For a finite, isotropic bath there is a crossing of the 
lowest energy levels at A c = hN/(N — 1). This is 
shown in Fig. [2fa), where we have plotted the energies 
E = {Hlmg) of the lowest levels for 7 = 1. The GS is the 
separable, spin-polarized state \N/2,N/2) for < A < A c 
(normal phase), and an entangled state \N/2, [hN/2X]) 
for A > A c (broken phase), where the function [x] rounds 
x to the nearest integer. For both phases the GS is 
non-degenerate for any bath size. Note that the criti- 



cal point depends on h, thus the bath can be directed to 
either phase by setting the magnetic field strength cor- 
respondingly greater or less than the intra-bath coupling 
strength. 

The GS of an anisotropic bath is a superposition over 
all Dicke states in the maximal spin subspace Sjy/2, but 
tends to the spin-polarized state \N/2, N/2) in the limit 
A — >• 0. The GS is non-degenerate for the normal phase, 
but becomes doubly degenerate in the thermodynami- 
cal limit for the broken phase. This is demonstrated in 
Fig. IHb), which shows that the energies of the lowest 
levels are almost identical for an anisotropic bath of just 
N = 30 spins. In the thermodynamical limit, the GS 
magnetization (per spin) in the z-direction is given by 

1 to \ f 1/2 for 0<A<ft ... 
Kgs = N {Sz) = { h/2X for \>h > < 4 > 

for all 7. See Figs. [2jc) and[2jd) for the magnetization of 
a large but finite bath with 7 = 1 and 7 = respectively. 

2. Analytical calculations with a partitioned LMG bath 

The high symmetry of the LMG bath allows us to eas- 
ily compute exact solutions of the system dynamics for a 
large number of spins, since the GS can be expressed in 
a reduced subspace S5 scaling at most as N + 1. How- 
ever, this representation of the GS is limited in that we 
can only consider operators that are symmetric across 
the entire bath ensemble. Fortunately, the LMG model 
provides some flexibility since we can make partitions of 
the bath, as shown in Fig. [1] without significantly in- 
creasing the scaling with system size, and thus can also 
consider non-symmetric operators. As an example, con- 
sider an operator that acts on, and is symmetric across, 
a total of k bath spins. We make a bipartite split of the 
bath into subsystems A and B with k and N — k spins 
respectively, and transform the original Dicke states as 

\(sa,sb)S, M) = ^2 \sA,m A ) <8> \s B ,m B ) 

niA,m B 

x (sA,m,A\ <8> {s B ,m B \S,M) , (5) 

where |saj to a) and \s B ,m B ) are Dicke states for sub- 
systems A and B; s A < k/2; and s B < (N — k)/2. Wc 
denote the tensor product space by Ss A <S)S>^~ k , where the 
superscripts indicate the corresponding number of spins. 
On the LHS of Eq. ([5]) we specify sa and s B because, in 
general, there are a range of allowed values these quan- 
tum numbers can take when coupling to form \S, M). 
However, for the case where S is maximal, which is our 
only concern, we have sa — k/2, s B — (N — k)/2 and 
sa + s B = S, reducing the Clebsch-Gordan coefficients 
in Eq. ((5|) from the general form in Ref. [43] to 
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/I/O Ifl/fir M/o ITU/O j|f\ X F m +k2 (k)F m +/yN _ k)2 {N - fc) 

{k/2,mA\® ({N -k)/2,mB\N/2,M) = 6 mA+mB , M \ ^ tttt , (6) 

y f m+n/2 \jy ) 

where F^a) = a\/[b\(a — b)l) are binomial coefficients. From Eqs. ([5]) and © we can construct a transformation 
operator that maps the subspace §n/2 to the tensor product space §JjS 2 ® given by 



JV/2 fc/2 {N-h)/2 

Tw/2^k/2,(N-k)/2 = ^ ^ <W+«Jb,-M^ 

M=-JV/2 m A = -fc/2 mB=-(N-k)/2 



< F mA+k/2 (k)F mB+{N _ k)/2 (N - fc) 



F M +N/2(N) 

x |fc/2,mA)® |(iV-fe)/2,m B )(iV/2,M| . (7) 



We apply this transformation operator to the representa- 
tions in §/v/2 of both the GS and -Hla/g, allowing us to 
consider operators that are symmetric across just fc spins 
without needing to revert to the full Hilbert space 1 



iigiJV 
'2 



we vary fc from 1 to N in order to consider different levels 
of exposure of the qubit to the bath. 



D. Characterizing decoherence 



B. Ising interaction 

The first interaction term we consider is of the Ising 
type 



Hi 



Vk 



k 

£< 

i=i 



(8) 



where e is the interaction strength. This commutes with 
the free qubit Hamiltonian (i.e. [Hj,Hq] — 0) and thus 
there will be no energy exchange between the qubit and 
the bath; the evolution will affect the off-diagonal terms 
of the state of Q only. The term above describes a situa- 
tion involving a total of fc out of N bath spins interacting 
with the qubit. By using the transformation operator in 
Eq. ([7]) we can vary the value of fc to study any scenario 
from a single coupled bath spin (fc = 1) to a completely 
coupled bath (i.e. the central spin model, fc = N), allow- 
ing us to evaluate how the evolution of the qubit changes 
as it is exposed to varying numbers of bath spins. We 
have introduced the factor 1/y/k so as not to artificially 
increase the interaction strength as we increase the num- 
ber of links, effectively normalizing the interaction term. 



C. LMG interaction 

With the second interaction we allow energy exchange 
between the qubit and the bath and thus introduce dissi- 
pation into the subsystem Q. A natural way of doing this 
is to consider the qubit interacting with the bath spins 
in the same way the bath spins do with each other, i.e. 



Hj = ±= (<T* <8> of + 7?" ® crV) . (9) 

v i—l 

The anisotropy parameter 7 used here is the same as 
that in the LMG Hamiltonian in Eq. ([3]). Once again, 



The main aim of this work is to characterize how the 
coupling to an LMG spin bath affects the decoherent evo- 
lution of the qubit Q. In the following subsections we 
introduce two measures that allow us to do this; namely, 
the average purity, and whether the induced decoherence 
is 'entanglement breaking'. Their calculation requires the 
use of an ancilla qubit and exploits the Jamiolkowski iso- 
morphism [44j], which is detailed in appendix [X] Impor- 
tantly, all the information pertaining to the decoherent 
evolution of the single qubit is mapped via the Jami- 
olkowski isomorphism to a two-qubit state p A for Q and 
the ancilla. From p A we can ascertain if and when the 
evolution is entanglement breaking, and also deduce the 
Kraus operators for the evolution to obtain the average 
purity. 



1. Average purity 

We use the purity of the reduced quantum state p(t) 
of the qubit, defined as P[p] = Tr(p 2 ), as a measure of 
its decoherent evolution in the presence of the bath. The 
qubit state at time t will be pure iff P(t) — 1, and mixed 
for P(t) < 1. We initialize the qubit in a pure state at 
t = so that any mixing of the state, i.e. decoherence, 
can be observed in its subsequent evolution described by 
the unital map e{p(0)}. In general, the purity depends 
on the exact details of the state at t = and thus we 
average over all possible initial pure states \tp) to obtain 
the average purity. This can be achieved once the Kraus 
operators Ai for the evolution are known via 45] 



P = Tr[ £ {|^)(V|} 2 



d(d + l) 



(10) 



where d — 2 is the dimension of the quantum subsystem. 
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2. Entanglement breaking 

A useful property to determine environmentally in- 
duced decoherence is whether it destroys any entangle- 
ment the primary subsystem Q has with an external sub- 
system. The superoperator A (see appendix |A"]) describ- 
ing the effect of the interaction is 'entanglement break- 
ing' (4(| if, when acting on a subsystem 6 of dimension db, 
the final state p{™ = 1 ® A a (p™) is separable for every 
(possible entangled) initial state p™ of the composite 
system of b and an external subsystem a of dimension 
d a (not necessarily a copy of b). Importantly, when the 
state p A obtained as described in appendix [A] is separa- 
ble, the superoperator A will be entanglement breaking 
since this implies the corresponding quantum operation e 
has a Kraus representation composed entirely of projec- 
tors. Separability can be determined , using the positive 
partial transposition (PPT) criterion [4^ [48], by comput- 
ing the minimum eigenvalue p(t) of the matrix (p A ) Ta ■ 
The induced decoherence will be entanglement breaking 
iff>(t)>0. 

III. RESULTS: ISING INTERACTION 

In this section we present results for the single spin 
decoherence induced by the Ising interaction with the 
LMG bath spins. Before doing so, it is useful to discuss 
exactly how decoherence arises in this system. 

Consider the qubit initially prepared in a generic su- 
perposition of its ground and excited states, such that 
the global state at time t = is 

|*(0)) = (o |0) + oi|l))® \GS) , (11) 

where the coefficients ao and a\ satisfy |ao| 2 + |&i| 2 = 1. 
The subsequent evolution under Ht for times t > is 
described by 

|*(t)) = aoe- lut \0) ® e-^+^GS) 

+ aie iut |l) <8> e- iH - l \GS) , (12) 

where S k z = Ei=i°fA and H ± = H R ± 2eS^/Vk 
are perturbed bath Hamiltonians. In general, the bath 
evolves to different states dependent on whether the 
qubit was initially in the state |0) or |1), and therefore 
Eq. (fT2|) is an entangled state of the two subsystems. 
After tracing out the bath we obtain the reduced density 
matrix for the qubit (in the eigenbasis {|0), |1)}) 

p(t) = Tr R {Mt))m)\}=( l p ^ ) (13) 
where 

Pi 2 (t)=a ate- 2 ^*(t) . (14) 

The relative populations of the ground and excited states 
of the qubit are conserved, but dephasing noise is intro- 
duced through the decoherence factor 

D(t) = {GSy^e-^-^GS) . (15) 



Note that the induced decoherence is independent of u. 
The purity of the state in Eq. (TIB")) is 

P[p(t)] = l-2\a Q \ 2 \a 1 \ 2 [l-L(t)} , (16) 

where L(t) = \D(t)\ 2 is a real quantity taking the form 
of a Loschmidt Echo (LE) [49| for the bath. Unlike the 
purity, L(t) is independent of the initial state of the qubit, 
aside from the fact that it is undefined for ao = 0,1. 
Thus , the average purity is trivially related to the LE via 
P(t) = [2 + L(t)}/3. 

By the term Loschmidt Echo we mean that L(t) quan- 
tifies the distance between the two bath states that 
have evolved from the GS under the two perturbed bath 
Hamiltonians H±. When the distance between these 
states is small, the entanglement between the qubit and 
the bath is low and the average purity is close to one. 
Conversely, when the LE is zero the two time-evolved 
bath states are orthogonal and the qubit is maximally 
entangled with the bath. In such cases, the qubit un- 
dergoes total dephasing to a completely mixed state and 
the average purity takes its minimum value of 2/3. Also, 
this type of dephasing noise is entanglement breaking iff 
I(i) = O0. 

A useful quantity to consider that will give insight into 
how the decoherence varies with link number k is the 
Fourier Transform of the decoherence factor in Eq. f| 15[) , 
given by 

X {E)=Y J (GS\<pi)(4>t\^)(^\GS)5{E-EJ +E+) , 

i,j 

(17) 

where ) are eigenvectors of the perturbed bath Hamil- 
tonians H±, and Ef 1 their eigenenergies. Contributions 
to x{E) occur at energy differences EJ — Ef correspond- 
ing to eigenstates \4>J) an d I'frt) that simulaneously con- 
nect to (i.e. are not distant from) each other and the 
bath GS. 



A. Single link, k = 1 

The simplest case to consider first is that of a single 
link between the qubit and the bath. In particular, for 
an isotropic bath (7 = 1) the GS in the normal phase 
(A < h) is the spin-polarized state \N/2,N,2), which is 
an eigenstate of both perturbed bath Hamiltonians H± . 
The bath remains in the GS for the entirety of the evolu- 
tion and simply induces an additional z-rotation on the 
initial qubit state through an angle 2et, thus there is zero 
dephasing. 

In the broken phase (A > h) we observe a loss of qubit 
coherence and find that the average purity has an oscil- 
latory behavior with time, as shown in Fig. [3Ja) . For an 
isotropic bath, we show in appendix iBlthat the GS in the 
broken phase decomposes into two terms when we move 
to the bipartite subspace §J , % ® ^N-iy 2 - ^he evolution 
of the bath under H± is further confined to a subspace 
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comprising only the two states of the decomposition, al- 
lowing us to describe the bath propagators e~ lH±t by 2 
x 2 matrices (see Eqs. (|B9j) and (|B10[1 ). The resulting 
expression for the purity is too complicated to be given 
explicitly however, we can extract the salient behavior 
from the form of the two propagators. The off-diagonal 
matrix elements of e -1 ^-* vary sinusoidally such that it 
periodically becomes equal to the identity matrix with a 
frequency r\- — 2yJ\ 2 + e(e + 2h) that is independent of 
bath size N. Similarly, e~__ varies periodically with 
a frequency t] + = 2^J\ 2 + e(e — 2h). When both prop- 
agators are the identity at the same time, the LE be- 
comes unity and there is a revival of full qubit coher- 
ence. This 'coherence time' t c must simultaneously sat- 
isfy rj-T c = 2Z_7T and rj + t c — 2l + n, where Z_ and are 
integers. In general, we observe fast oscillations of the 
average purity determined by the higher frequency T)-, 
which are contained within an envelope function that is 
periodic with a frequency rj+. High average purities of 
P(t, A) > 0.98 are obtained at times t = W r (I an in- 
teger), where r r = Tr/i] + is the approximate 'rephasing 
time' corresponding to the first maximum of the envelope 
function. This is plotted in Fig. |3£b) for various values 
of the interaction strength e. In the limit e <C h, the 
frequencies r}± are approximately equal and a full revival 
of qubit coherence occurs at a time r c ~ t t = tt/A. For 
e = h the rephasing time is asymptotic at A = h, and 
thus its sensitivity to small changes in the magnetic field 
allows for an observation of criticality. However, critical- 
ity should be avoided if one wishes to use the system for 
storage of the qubit state so that rephasing can occur on 
short-time scales. 

The behavior described above for the coupling to an 
isotropic bath is independent of bath size N. Thus, the 
revivals of qubit coherence are not a finite size effect and 
occur in the thermodynamical limit. This independence 
with N arises because the bath Hamiltonian takes a par- 
ticularly simple form for 7=1 that is diagonal in the 
eigenbasis spanned by Dicke states \S,M), ultimately 
restricting the evolution to just a two-dimensional sub- 
space. As a result, the decoherence factor has a maxi- 
mum of four Fourier components, as shown in Fig. EJc); 
there is only a single component in the normal phase 
where no dcphasing occurs. Interestingly, when we con- 
sider anisotropic baths such that the evolution utilizes 
the entire §i/2®§(iv^i)/2 subspace, we find there are still 
only n — 4 significant Fourier components to the deco- 
herence factor. Significant Fourier components are those 
of greatest amplitude that reproduce 98% of the signal. 
Also, we observe a rapid convergence of the average pu- 
rity to the isotropic bath result with increasing bath size. 
In Fig. EJd) we have plotted the average purity across 
criticality at time t = 5/h for various bath anisotropics 
and a bath of AT = 100 spins. Away from criticality, 
the average purity for all bath anisotropics is within 5% 
of the isotropic result for a bath of this size. Close to 
criticality the convergence is slower, and we observe os- 
cillations in the average purity in the normal phase that 
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FIG. 3. Ising interaction, single link (k — 1): (a) the av- 
erage purity plotted against A and t. (b) The approximate 
rephasing time r r plotted against A for various values of e. (c) 
The Fourier Transform of the decoherence factor as a func- 
tion of energy and A. We have used Gaussian envelopes at 
each Fourier component for ease of viewing, (d) The average 
purity plotted against A, at t — 5/h, for various values of 7. 



disappear in the thermodynamical limit. 



B. Completely connected, k — N 

We now consider the opposite extreme where the qubit 
is coupled to all spins in the bath, i.e. the central spin 
model. For an isotropic bath the average purity is triv- 
ially equal to one for all A; the interaction term commutes 
with the bath Hamiltonian and thus the qubit and the 
bath cannot become entangled. The qubit experiences a 
z-rotation through an angle eNt in addition to that of its 
free evolution [50(. 

Figure [UJa) shows a color map of the average purity 
against A and t for an anisotropic bath (7 = 0) with 
e = 0.25/i. At this low interaction strength the aver- 
age purity is qualitatively similar to the single link sce- 
nario; there is almost zero dephasing in the normal phase, 
whilst oscillations occur within an envelope function in 
the broken phase. The rephasing time t t is asymptotic 
at criticality but contrary to the single link scenario we 
find it is independent of e. As would be expected for this 
description the Fourier Transform of the decoherence fac- 
tor for low interaction strengths is qualitatively similar to 
that shown in Fig. [3^c); there is one Fourier component 
with a large amplitude accompanied by a few (n < 15) 
of smaller amplitude. 

As the interaction strength is increased above e = h the 
behavior becomes markedly different. The maxima in the 
average purity are suppressed and, as shown in Fig. Sfb) 
for e = 5h, high average purities are no longer achieved on 
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FIG. 4. Ising interaction, completely connected (k — N): (a) 
a colour map of the average purity plotted against A and t for 
e = 0.25/i. The map has a threshold at P = 0.99 to improve 
clarity, (b) The average purity plotted against A and t for 
e = 5h. (c) The Fourier transform of the decoherence factor 
as a function of energy and A for e = 5h; again, we have 
used Gaussian envelopes at each Fourier component, (d) The 
number of Fourier components n of greatest amplitude that 
reproduce 98% of the signal, plotted against A, and for various 
values of e. 



short-time scales, except in the limit A — > 0. In addition, 
away from the maxima the LE is close to zero indicating 
that the evolution is almost entanglement breaking. The 
Fourier Transform of the decoherence factor for e = 5h, 
shown in Fig. EJc), contains a large number n of signif- 
icant components but no single one of high amplitude, 
resulting in the suppression of the maxima. In Fig. SJd) 
we have plotted n across criticality for various values of e. 
The number of significant Fourier components increases 
rapidly as we move from the normal to the broken phase, 
and contrary to the single link scenario increases at all 
intra-bath coupling strengths with e. 

Interestingly, for any e the energy spread and number 
of components rapidly converge to finite values with in- 
creasing bath size. Thus, although at high interaction 
strengths there is a decay of qubit coherence on short- 
time scales, full revivals will eventually occur for any e 
even in the thermodynamical limit, with the coherence 
time r c independent of N for large baths, N > 100. Fi- 
nally, we note that the behavior described in this sub- 
section is only qualitatively the same for different bath 
anisotropics and does not tend to the isotropic behavior 
in the limit 7 — > 1. 



C. Multiple links 

In the previous two subsections we have seen approxi- 
mate rephasing occurs on short-time scales for any inter- 



action strength when the qubit is connected via a single 
link to the bath, whilst for a completely connected qubit 
rephasing is suppressed at high interaction strengths. We 
now examine how this transition occurs as we gradually 
increase the number of links from the scenario k = 1 to 
k = N . We consider an intra-bath coupling strength of 
A = 3h such that we are away from criticality and rephas- 
ing is possible on short-time scales for any link number. 

In Figs. 03a) and [5jc) we have plotted the maximum 
value of the average purity for rephasing within a time 
t < 10//i, as a function of the interaction strength e, and 
for a range of link numbers k. The bath size is N — 50. 
High average purities of P > 0.98 are achieved for any 
link number below e ~ Ah, whilst above this interaction 
strength a high average purity is only guaranteed for k < 
3. At higher link numbers we observe the suppression of 
maxima up to e ~ lOh as discussed in Sec. IIIIBI for a 
completely connected bath. Interestingly, above e = lOh 
the maxima in the average purity increase back towards 
one, with the rate of increase slowest for link numbers k ~ 
30-40. This behavior is prevalent when we consider the 
number of significant Fourier components n as a function 
of link number and interaction strength, shown in Fig. 
0Jb). For a few link numbers, we have n < 10 for all 
e allowing approximate rephasing to occur. For higher 
link numbers, the number of components increase with 
e up to a maximum value at e ~ lO/i, before decreasing 
once more at higher interaction strengths. The number of 
components also has a maximum at k ~ 35 and therefore 
decreases with large k, but albeit is still relatively large 
for a completely connected qubit. 

The restriction of the number of significant Fourier 
components can be explained by considering the form of 
the perturbed bath Hamiltonians H± — Hr ± 2eS^ / \fk. 
At low interaction strengths e < h, the perturbing terms 
are small and the eigenstates of the two Hamiltonians 
are similar. There are approximately (k + 1)(N — k + 1) 
terms (<pf\4>~) that are non-zero, but only a few of these 
will connect to the GS and thus x(-E) will have only a 
few components. As the interaction strength is increased 
above e = h, the two sets of eigenstates differ meaning 
that there are now (k + 1) 2 (N — k + l) 2 non-zero-terms 
((f>f\4>~), each of which will be of smaller magnitude. 
There is a larger set of pairs of eigenstates to connect 
to the GS leading to many more components. When the 
interaction strength is increased further such that the 
second terms dominate in H±, the eigenstates again be- 
come similar and the number of components asymptoti- 
cally decreases to one in the limit e — > 00 . The maximum 
in the number of components with link number arises be- 
cause the subspace §fe/ 2 ( 8 ) ^(jv r - fc fc)/2 ^ s iar g es t for k — N/2, 
however the maximum is skewed towards k — N so that 
a completely connected qubit still has a relatively large 
number of components compared to the case k = 1. 

Finally, we note that although the maxima in the aver- 
age purity depend on link number, the minima do not, as 
shown in figure [SJd). The minima decrease towards, but 
do not reach, the minimum value of 2/3 (corresponding 
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FIG. 5. Ising interaction: (a) the maximum value of the 
Loschmidt Echo at A = 3 for rephasing within a time t < 
10/h, as a function of e, and for a range of k. (b) The number 
of Fourier components n of greatest amplitude that reproduce 
98% of the signal as a function of e and k. (c) same as (a) but 
for only a few links, (d) The minimum value of the Loschmidt 
Echo within t < 10/h at A = 3, plotted against e. 

to an LE of zero) with increasing interaction strength and 
so the interaction is never formally entanglement break- 
ing. 




Ar=100,fc = jV,e = 0.25/i,7 = ' N = m,k=l,X = 3h 




FIG. 6. LMG interaction, single link (k — 1): (a) the average 
purity plotted against A and t for 7 = 1. (b) The rephasing 
time r r for the bath in the normal phase (A = 0.5), as a 
function of e and for various values of 7. (c) A color map of 
the average purity plotted against A and t for e = h. The 
map has a threshold at P = 0.75 to improve clarity, (d) The 
maximum value of fj,(t) within t < 10/h at A = 3h, as a 
function of e and for various values of 7. Values above zero 
indicate that the evolution is entanglement breaking. 



in the following results using the quantity ji(t) discussed 
in Sec. HTD21 



IV. RESULTS: LMG INTERACTION 

We now present the results for the qubit interacting 
with the bath spins via the dissipative LMG interaction. 
Contrary to the Ising interaction, the relative populations 
of the ground and excited states of the qubit can vary in 
time in addition to the coherence. This means the aver- 
age purity can take a minimum value of 1/2, but more 
significantly it is now dependent on the energy difference 
uj between the two levels of the qubit. In the following, 
we consider the scenario [5l| in which the qubit is subject 
to an identical external field as the bath spins, thereby 
setting lo — h. In tests where ui was varied decoherence 
was suppressed in all cases as uj was increased, thus the 
scenario considered is in the worst-case regime. 

One important consideration for the LMG interaction, 
in contrast to the Ising interaction, is that the induced de- 
coherence is non-trivially dependent on the initial qubit 
state. This means that the purity for a given initial state 
can have significantly different behavior to the average 
purity, and indeed can be unity when the average purity 
is not. It is interesting in the following results that we ob- 
serve the average purity oscillating away from and back 
to unity, indiciating universal rephasing at certain times 
for any given initial qubit state. The evolution is also en- 
tanglement breaking for certain periods, which is shown 



A. Single link, k = 1 

In Fig. Eta) we have plotted the average purity for 
the qubit connected via a single link to an isotropic bath 
(7 = 1) with an interaction strength e = h. For the 
LMG interaction the average purity is not equal to unity 
in the normal phase for an isotropic bath, but instead 
we observe approximate revivals with a rephasing time 
T r that is independent of intra-bath coupling strength 
A. The rephasing time increases rapidly with interaction 
strength and is almost independent of bath anisotropy, 
as shown in Fig. \$[b). The interaction is never entangle- 
ment breaking for any 7 or e in this phase. 

As for the Ising interaction, we observe approximate 
revivals in the broken phase for all 7 and e. The rephas- 
ing time r r behaves in a similar fashion to the Ising case 
with a single link: it is almost independent of 7; increases 
with interaction strength up to e = h where it is asymp- 
totic at criticality (see Fig. Etc)); and decreases at higher 
interaction strengths. In Fig. Etd) we have plotted the 
maximum value of fi(t) within a time t < 10/h for various 
bath anisotropies. We observe that the induced decoher- 
ence is periodically entanglement breaking in time for 
7 > 0.1 and interaction strengths e > 8h. 

For the LMG interaction we again find that the average 
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FIG. 7. LMG interaction, completely connected (k = N): 
(a) and (b), the average purity plotted against A and t for 
7 = and interaction strengths e = h and e = 5h respectively, 
(c) and (d), the maximum value of fi(t) within t < 10/h at 
A — 0.5h and A = 3h respectively, as a function of e and 
for various values of 7. Values above zero indicate that the 
evolution is entanglement breaking. 
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FIG. 8. LMG interaction: (a) and (b), the maximum value of 
the purity at A = 0.5h and A = 3/i respectively, for rephasing 
within a time t < 10//i, as a function of interaction strength 
e and for various values of k. (c) and (d), the maximum value 
of fi(t) within t < 10/h at A — 0.5h and A = 3h respectively, 
as a function of e and for various values of k. Values above 
zero indicate that the evolution is entanglement breaking. 



purity is independent of bath size for large N. Thus, not 
only would we observe revivals in the thermodynamical 
limit, but interestingly they will occur after periods when 
the induced decoherence is entanglement breaking. 



ment breaking on short-time scales for all 7 above e ~ 6h 
in this phase. However, at higher interaction strengths 
this is no longer the case; for 7 = the decoherence is 
not entanglement breaking for e > 23h (see Fig. IHJd)). 



B. Completely connected, k = N 



C. Multiple links 



Figures [TJa) and [T](b) show the average purity, as a 
function of intra-bath coupling strength and time, for a 
completely connected qubit with 7 = and interaction 
strengths e = h and e — 5h respectively. For this bath 
anistropy there is a decay of the average purity in the nor- 
mal phase even at low interaction strengths. Note that 
the decay also occurs at A = for which the bath spins 
do not interact with each other, and extends across criti- 
cality into the broken phase at high interaction strengths 
e > h. Rephasing only occurs on short-time scales in the 
limit 7 — > 1, with the rephasing time r r independent of 
bath size for large baths, N > 100. The interaction is 
now entanglement breaking in this phase on short-time 
scales except in the limit 7 — > 1, as shown in Fig. [TJc). 

Approximate revivals occur in the broken phase away 
from criticality with a frequency much greater than that 
of the corresponding single link scenarios. Contrary to 
the Ising case, the rephasing time r r increases with inter- 
action strength and converges less rapidly towards a finite 
value with increasing bath size. Therefore, the rephasing 
time can only be considered independent of bath size for 
very large baths, N > 1000. Figure EJd) shows that 
the induced decoherence is now periodically entangle- 



As was the case for the Ising interaction, for the dis- 
sipative LMG interaction we observe a change in the 
qubit 's behavior as we move from a single link to fully 
connected to the bath. In this subsection we discuss how 
this transition occurs as we vary the link number k for 
7 = 0. For this bath anisotropy approximate rephasing 
occurs on short-time scales in both phases for a single 
link, but is suppressed in the normal phase for k = N . 
Also, the decoherence is never entanglement breaking for 
k = 1 with 7 = 0, but is at high interaction strengths 
for k — N. We consider intra-bath coupling strengths 
away from criticality of A = 0.5/i for the normal phase 
and A = 3/i for the broken phase. 

Figures [5]Ja) and[5Jb) show the maximum value of the 
average purity for rephasing within t < 10/ h with a bath 
of N — 50 spins. Both figures are qualitatively similar to 
Fig. OJa) for the Ising interaction; approximate rephas- 
ing occurs for a single link at any interaction strength, 
whilst the maxima in the average purity are suppressed 
at high interaction strengths for greater link numbers. 
In the broken phase, high average purities of P > 0.98 
are achieved for any link number below e ^ Ah, whilst 
in the normal phase this isn't the case and there is an 
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immediate decay of the average purity above e = for 
k ^ 1. As for the Ising interaction, the suppression of 
maxima is reduced as the interaction strength is further 
increased, with the rate of increase slowest for link num- 
bers k ~ 30-40. This suggests that our analysis discussed 
in Sec. IIII C\ where the induced decoherence was limited 
by the size of the subspace required to describe the par- 
titioned bath states, has a broader applicability than the 
Ising interaction alone. Rephasing depends on the num- 
ber of significant Fourier components for, in this case, 
the average purity. This quantity is reduced in the lim- 
its e <C h and e 3> h for which the interaction term is 
a small perturbation on the bath Hamiltonian and vice 
versa respectively. 

In Figs. [8jc) and|Sfd) we have plotted the maximum 
value of p(t) within t < 10 /h for the normal and broken 
phases respectively. In the normal phase entanglement 
breaking occurs for k ^ 1 at interaction strengths e > 
0.8h, whilst in the broken phase the induced decoherence 
is entanglement breaking for k ^ 1 within the range h < 
e < 23h. It is only the case of coupling via a single 
link that the qubit is immune to entanglement breaking 
within the period t < 10/ h. 



V. CONCLUSIONS 

To summarize, we have investigated the quantum evo- 
lution of a single qubit coupled to a Lipkin-Meshkov- 
Glick bath as a model for decoherence in solid state quan- 
tum memories. The bath, which exhibited a second-order 
quantum phase transition, was highly symmetric and al- 
lowed for exact calculations of the system dynamics for 
large system sizes, N ~ 100 spins. Further, partitioning 
of the bath was possible without significantly increasing 
the computational complexity, which allowed us to de- 
termine the effect of increasing exposure of the qubit to 
the bath spins. Decoherence of the qubit was quantified 
using the average purity and determining if and when the 
evolution destroys any entanglement the qubit may have 
with an external subsystem. 

For the qubit interacting via an Ising interaction with 
just a few bath spins, we observed zero dephasing in the 
normal phase and almost complete revivals of qubit co- 
herence in the broken phase. The rephasing time, which 
was sensitive to criticality in the bath, was independent 
of bath size for large baths and thus faithful qubit storage 
would be possible in the thermodynamical limit. In gen- 
eral, as the number of links was increased the revivals in 
the purity were suppressed, which was discussed in terms 
of the number of components in the Fourier Transform of 
the decoherence factor. Decoherence was suppressed by 
the restrictive size of the subspace required to describe 
the partitioned bath states, which was smallest for just 
a few links. Further, it is likely that revivals were ob- 
served in certain scenarios, even in the thermodynamical 
limit, because the symmetry of the bath constrained the 
dynamics to a subspace whose dimension grows only lin- 



early with system size. Such revivals could possibly oc- 
cur in other systems possessing exchange symmetry that 
are not as easily amenable to analysis as the LMG model. 
However, they may be difficult to observe experimentally 
due to a potential lack of the required bath symmetry in 
real samples. 

In contrast to the Ising interaction, for the dissipative 
XY-type (LMG) interaction between the qubit and the 
bath, decoherence was dependent on the energy differ- 
ence between the two levels of the qubit and additionally 
was observed in both bath phases. Revivals were found to 
occur for certain parameter regimes and just a few links 
between the qubit and the bath. These were once again 
suppressed as the link number was increased, broadening 
our analytical discussion for the case of the Ising interac- 
tion in terms of the restrictive size of the bath's Hilbert 
space. Interestingly, the revivals occurred after periods 
when the evolution was entanglement breaking, indicat- 
ing a constant transfer of quantum information back and 
forth between the qubit and bath. Our results were based 
on a worst-case regime in which the qubit was subjected 
to an identical transverse field to the bath and in general 
was found to be less robust against decoherence than the 
Ising interaction. Thus, systems with Ising interactions 
coupling only to a small neighborhood of the environment 
would perform better as quantum memories. 
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Appendix A: Quantum Operations Formalism and 
the Jamiolkowski Isomorphism 

A powerful tool for considering the evolution of a quan- 
tum system is the quantum operations formalism (l|. If 
the initial state of the quantum system is described in 
terms of a density operator p, the subsequent evolution 
causes a transformation to a final state given by the map- 
ping p — > e{p}. The linear, completely positive map e 
is known as a quantum operation. For a closed system 
that implements a particular unitary operation U, the 
quantum operation is simply e(p) = UpU^. In the same 
context, the evolution of an open quantum subsystem 
interacting with some environment R, e.g. a single spin- 
1/2 interacting with an LMG spin bath, can be described 
as 

E(p) = Tr R [U T (p®p R )Ul] , (Al) 

where we consider the total system as closed and obtain 
the final reduced density operator for the principal sub- 
system by tracing out the state of the environment (note 
that we assume an initial product state for the subsystem 
and environment). 
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Although an elegant description, calculating the final 
state in the above way is often computationally difficult 
due to the large size of the total Hilbert space. We can 
instead express Eq. (|A1[) explicitly in terms of the prin- 
cipal subsystem's Hilbert space H by 



e{p)=Y j A iP A\ 



(A2) 



where the A%, which act on H, are known as Kraus op- 
erators and d is the dimension of %. Once the Kraus op- 
erators are known for a particular type of evolution, one 
can easily obtain the final state of the open subsystem 
given any initial state. Importantly, the Kraus operators 
reveal the nature of the noise induced by the coupling 
to the environment. For a unital map, which is always 
the case for the noise described by the coupling to an 
environment, the Kraus operators satisfy AiA\ = 1. 

An equivalent description of the evolution in Eq. (IA2() 
can be made in terms of a superoperator A. To explain 
this more fully, consider the Hilbert space H spanned 
by basis states {\i) | i = 0, • • ■ , d — 1}. We can expand 
any density operator p for the system in the operator 
basis | i,j = 0, - ■ ■ ,d— 1}, with its corresponding 

matrix elements pij contained in a c? 2 -dimensional vector. 
The superoperator A can then be described by a d 2 x d 2 
super-matrix with elements 



are related to those of A by dpf k ^ — h-ijki , where 

P A = E Piw\i)U\®\k)(l\ • (A6) 

ijkl=0 

Thus, using the Jamiolkowski isomorphism, it is sufficient 
to calculate p A to determine the superoperator A, and 
then trivially the Kraus operators Ai, for any quantum 
operation e(p). 
Appendix B: The Loschmidt Echo for a single link 
and an isotropic bath 

The GS of an isotropic, ferromagnetically coupled bath 
is the Dicke state \N/2,M), where M = N/2, [hN/2X] 
in the normal and broken phases respectively. For Af>l 
the latter can be approximated as M = hN/2X. When 
we partition the bath into a single and N — 1 spins the 
Dicke state \N/2,M) decomposes into two terms 



I JV 
I 2 ' ' 



,M> = Cl/2 ||,i>®|f -i M-i> 

+ c_i/ 2 |i-£)® |f M + i) , 
where the Clebsch-Gordon coefficients are given by 



(Bl) 



- M 

2 ' 



1 I JV 

2 I 2 



C-l/2 - (|, 



(f 



1 I N 
2 



f,M) = 



M 
N 



M 
N 



(B2) 



d-i 



(A3) 



fei=0 



Note that A can be inferred from the Kraus operators 
and vice versa via 



A = ^ A * ® A i 



(A4) 



The Jamiolkowski isomorphism |44j | exploits an initial 
setup of two copies of the principal subsystem, such that 
the subsequent evolution transfers all the information 
about A to the final quantum state. The two copies 
a and b are prepared in the maximally entangled state 
|* + ) = J2i=o N) ® N) an d th° superoperator A is ap- 
plied to b as 



d-i 

(l®A){|¥ + ><*+|} = -£ ® A\i)(j\ = P A . (A5) 

ij=0 



The matrix elements of the resulting density operator p 



Evolution of the GS under the perturbed bath Hamilto- 
nian H_ = Hr — eaf for this scenario is restricted to the 
two states on the RHS of Eq. (|Blj) . Using these states 
as a basis, the action of H- on the Dicke state | N/2, M) 
can be summarised by the matrix equation 



H- |f ,M) = h 
where 



e/h 



P 



P a + C + e/h 



Cl/2 
C-l/2 



A 



2hN 



{(N - l) 2 - (2M - l) 2 } - 2M 



T\ l 



AM 2 



4AM 
hN 



For N 3> 1, the latter two quantities become 



C = 2. 



(B3) 

(B4) 

(B5) 
(B6) 



(B7) 
(B8) 



By firstly diagonalizing the matrix representation of H— 
in Eq. (|B3|) . we can find a similar representation for 



,-iJ?_t 



in this basis, given by 
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-m.t _ e-^+W* ft,- cos (q-t/2) + 2i(h + e) sin -2i/3/ l sin(r ? _t/2) , 

77_ V -2i/3/isin(r?_*/2) r?_ cos (»7-*/2) - 2i(/i + e) sin {r)-t/2)} ' 1 



where ?7_ = 2-^/A 2 + e(e + 2/i). Similarly, for e tH+t we obtain 

^or+t _ e-^+W f v+ cos (r? + t/2) + 2i(h - e) sin (r? + t/2) -2H3hsm(r, + t/2) * 

jy + I -2i/3/isin(r?+V 2 ) V+ cos 0?+V 2 ) - 2i(/i - e) sin (r?+i/2W ' 



where ry + = 2-\/A 2 + e(e — 2/i). The LE for the broken 
phase (A > h) can now be calculated by multiplying the 
hermitian conjugate of Eq. (|B10|) by Eq. (|B9[) and taking 



the expectation value with respect to the GS, i.e. 



L(t) 



(Ci/2,C_i/ 2 



)e^ 



iH-t 



Cl/2 
C-l/2 



(Bll) 



The resulting expression for the LE is complicated and 
need not be written here explicitly. Importantly, we can 
see revivals of full qubit coherence occur at times t = It c 
that satisfy both ?7_t c = 2Z_7r and ??+t c = 2Z + 7r, where 
I- and 1+ are all integers; at these times the matrices in 
Eq. (|B9f and Eq. (|B10j) are both equal to the identity. 
In the limit e <C h the two frequencies are approximately 
equal and the coherence time is r c = ir/X. 
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